Suppression of 1/ f a noise in one-qubit systems 
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We investigate the generation of quantum operations for one-qubit systems under classical noise 
with l/f a power spectrum, where 2 > a > 0. We present an efficient way to approximate the 
noise with a discrete multi-state Markovian fluctuator. With this method, the average temporal 
evolution of the qubit density matrix under 1/ f a noise can be feasibly determined from recently 
derived deterministic master equations. We obtain qubit operations such as quantum memory and 
the NOT gate to high fidelity by a gradient based optimization algorithm. For the NOT gate, the 
computed fidelities are qualitatively similar to those obtained earlier for random telegraph noise. 
In the case of quantum memory however, we observe a nonmonotonic dependency of the fidelity on 
the operation time, yielding a natural access rate of the memory. 



I. INTRODUCTION 



In solid-state realization of qubits, material specific 
fluctuations typically induce the major contribution to 
the intrinsic noise. Much effort has been focused on 
the preservation of the state in a quantum memory in 
the presence of l/f a noise since this is a ubiquitous 
form of noise encountered in solid-state qubit applica- 
tions [l|, 0, Q. Both charge and spin qubits are sus- 
ceptible to noise of this form. For Josephson junctions, 
both charge noise [3, 0] and critical current noise 0, @] 
have been measured to have l// a power spectral densi- 
ties. Similar charge fluctuations are responsible for the 
well-known 1/ f a nature of low frequency noise in sin- 
gle electron transistors Background charge fluctua- 
tions resulting in 1// Q noise spectra are considered to 
be the most important source of dephasing in Joseph- 
son junction qubits 0, H, 0|. Spin qubits such as those 
formed from donor spins in semiconductors are suscepti- 
ble to nuclear spin noise deriving from dipolar coupling 
between environmental nuclear spins. The nuclear spin 
bath couples to the donor spins by hyperfine interac- 
tions, which renders the dynamics of the nuclear spins 
to cause dephasing. Recent calculations for a phospho- 
rus donor in silicon show that the high frequency com- 
ponent of the nuclear spin noise is approximately de- 
scribed by a l/f a power spectrum [lfj]. Electron spin 
qubits implanted into silicon are also affected by 
relaxation of dangling bonds deriving from oxygen va- 
cancies at the Si/Si02 interface. This gives rise to a 
magnetic noise with a 1//° spectrum that is the dom- 
inant mechanism for phase fluctuations of donor spins 
near the surface Another form of noise closely re- 

lated to 1// Q noise is random telegraph noise (RTN), 
which arises from coupling of individual bistable fluctu- 
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ators to a qubit @, O, B [H, [H, [T3, [H, [T§| . 

Several approaches to suppress decoherence based on 
pulse design have been proposed in the literature. Among 
them, dynamical decoupling schemes average out the un- 
wanted effects of the environmental interaction through 
the application of suitable control pulses [13, HH . Appli- 
cation of these schemes often involves hard pulses with 
instantaneous switchings and unbounded control ampli- 
tudes, resulting in a range of validity restricted to time 
scales for which the pulse duration is much less than the 
noise correlation time [H, US] . In Ref. [13] , a direct pulse 
optimization method restricted to bounded control pulses 
was developed for implementing one-qubit operations in a 
noisy environment. This initial work on noise suppression 
addressed the example of a single qubit system under the 
influence of classically modeled random telegraph noise, 
such as might arise from a single bistable fluctuator. 

In this paper, we extend the work of Ref. [13] to 
the physically relevant situation of l/f a noise where 
2 > a > 0. This kind of noise is known to result jfor ex- 
ample, from a set of bistable fluctuators [25l. I26L I27L |28| . 
i.e., RTN sources. We investigate two ways to approxi- 
mate the l/f a noise for computer simulations, namely, 
the sum of independent RTN fluctuators and a single 
discrete multi-state Markovian noise source. We show 
that the single fluctuator provides a much more efficient 
way to model l// a noise than independent RTN fluctu- 
ators. Furthermore, the average temporal evolution of 
the density matrix under this Markovian noise can be 
exactly described by a set of deterministic master equa- 
tions derived in Ref. (HI- Using this approach, we avoid 
the heavy computational task arising from the numerical 
evaluation of the density matrix averaged over a large 
number of different sample paths of the noise as com- 
puted in Ref. [13] • This framework will not only signif- 
icantly accelerate the convergence of the control pulse 
sequence optimization, but also allows further theoreti- 
cal analysis. Using these master equations, we employ 
gradient based optimization procedures to obtain pulse 
sequences that suppress 1//" noise for quantum mem- 
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ory and for a NOT gate. Comparisons with composite 
pulses designed to eliminate systematic errors and with 
refocusing pulses demonstrate that the numerically opti- 
mized pulse sequences yield the highest fidelities. 

The remainder of this paper is organized as follows. 
In Sec. HU we show how to efficiently approximate the 
1//" noise by a multi-state Markovian fluctuator. In 
Sec. IIIIl we define the fidelity of qubit operations, re- 
view the master equations describing the average evolu- 
tion of the qubit density matrix in the presence of the 
noise and describe the numerial optimization procedure. 
Sections |IV] and [V] present optimized control pulse se- 
quences and the achieved fidelities for quantum memory 
and for the NOT gate, respectively. Finally, Sec. lVII con- 
cludes and indicates further applications of the method. 



II. ONE-QUBIT SYSTEM SUBJECT TO l/f a 
NOISE 

We consider a one-qubit system described by the effec- 
tive Hamiltonian 



H = ^a(t)a x + ^T)(t)a z , 



(1) 



where alt) 6 [— a max , a max ] is the external control field 
applied along the x direction and r}{t) is the classical 
noise signal perturbing the system along the z direction. 
The noise source rj{t) can be characterized by its auto- 
correlation function 

C(t) = (v(O)v(t)) = Jim i / 1 v(s)v(s + f) ds, (2) 



-T/2 



the Fourier transformation of which defines the noise 
power spectral density as 



S(f) 



C(t)e- i27rft dt. 



(3) 



For a single RTN source with the amplitude A and cor- 
relation time r c , the autocorrelation function is given 
by 



C RTN (t) = A 2 e - 2 l t l/-, 
and the corresponding power spectral density by 



<Srtn(/) — 



1 + (7r/r c ) 2 



(4) 



(5) 



A standard way to simulate 1// Q noise is to use 
an ensemble of K independent uncorrelated RTN pro- 
cesses P, HH, Let r]k(t) be a symmetric RTN signal 
switching between values — and Afc with the corre- 
lation time r k = 1/7^, where jk is the transition rate 
between the two states. The total noise process appears 
in the Hamiltonian ([!]) as rj(t) = 2fe=i Vk{t)- Since the 



RTN sources are independent, Eqs. (J2j) and {4]) yield the 
autocorrelation function 



K 



K 



k=l k=l 

and the corresponding power spectral density is given by 



(7) 



Introducing the density of transition rates 3(7) and ex- 
pressing the noise strength A as a function of the tran- 
sition rate, we can replace the summation in Eq. J7]) by 
an integration, which yields 



S(f) 



7m " A 2 ( 7 )g( 7 )7 
m ,„ 7 2 + W) 2 



d 7 > 



(8) 



where 7 m j n and 7 ma x are minimal and maximal transition 
rates, respectively. Provided that 



A 2 ( 7 ) ff ( 7 ) = 2A/ 7 , 



(9) 



where A is a constant, the power spectral density in 
Eq. ((HI) becomes [H] 



A 



arctan 



In 



7n 



7T/ 

< Iff < In 



arctan 



In 
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Thus Eq. (JTUJ) yields an approximation to the 1// power 
spectrum. To generate a general l/f a power spectral 
density for 2 > a > 0, we can choose 



A 2 ( 7 ) 5 ( 7 ) = 2Ar a 



(11) 



as shown in [271 ] . 

Although the above method yields a valid approxima- 
tion for the l/f a spectrum, it is computationally ineffi- 
cient. In particular, the number of distinct noise states 
increases exponentially with the number of RTN fluctua- 
tors K, i.e., the number of terms in the sum of Eq. ([7]) ap- 
proximating the l// a noise. Since the size of the differen- 
tial equation system describing the average qubit dynam- 
ics increases linearly with the number of noise states [29| . 
in practice one has to restrict the computation to a rather 
small number of independent RTN fluctuators. 

To overcome this problem, we present a conceptually 
different way of generating the desired l/f a noise spec- 
trum using a single multi-state Markovian fluctuator. 
Consider a continuous-time Markovian noise process with 
M discrete noise states. Let T^j denote the transition 
rate from the jth state to the fcth one. In order to pre- 
serve total probability, we must have 



M 



for all k = 1,2,. 



,M. 



(12) 



Let us assume that the transition rates are symmetric, 
i.e., T = T T . Under this assumption the noise process 
has a steady-state solution in which the different noise 
states are equally probable. In order for the noise to be 
unbiased, i.e., (77) — 0, the amplitudes bk associated with 
the noise states must satisfy 



M 



b » = o. 

fe=i 

Thus the autocorrelation is given by 

C(t) = ( V (t) V (0)) = ^b T e r \%. 



(13) 



(14) 



Since T is symmetric, we can diagonalize it with an or- 
thogonal matrix V as T = VAV T , where the real diagonal 
matrix A = diag{Xk}kLi carries the eigenvalues of T in 
a descending order. Defining \ '■— ~Ju^ T ^' we rewr ^ e 
Eq. (JT4J) in the form of Eq. © as 



C(t) = X 



M 

''• A|i| x = E^ Afc|t| . 
fe=i 



(15) 



In order to use this multi-state Markovian fiuctuator to 
approximate 1// Q noise, we have to choose the eigen- 
values Afe and the amplitudes Xfc such that Eq. I|lip is 
fulfilled. Moreover, we must construct the orthogonal 
matrix V such that T = VAV T satisfies Eq. ([12]) , the 
amplitudes bk satisfy Eq. (fl3|l . and the off-diagonal ele- 
ments of r must be non-negative. 

One way to satisfy these requirements is to pick an 
integer m > 2 and set M — 2 m and to choose the eigen- 
values as 



{A fc }£l 



"2{0, 7 min , 7 min + S, 7 m in + 2(5, ... , 7max}, 



where 7 max = (M - 2)6 + 7 min and < S < 7 m i n - Hence, 
the distribution of the transition rates g(-f) is uniform 
on [7 m in, 7max]- Then we set V — H® m , where H is the 
Hadamard matrix 



Explicit calculation shows that these choices ensure that 
Eq. HI]) is satisfied. To fulfill Eqs. (IT ]) and (fl3 ]l . we set 

Xi = and Xk = Ik*^ for A; = 2, ... , M, where j k 
is equal to 7 m i n + (k — 2)8. It can be shown that this 
construction will also produce transition matrices V with 
non-negative off-diagonal elements. Hence we have pro- 
vided an efficient way to implement l/f a noise. Note 
that the M-state Markovian fiuctuator, Eq. (fl5l) . cor- 
responds formally to Eq. |(6|) with M — 1 non-vanishing 
RTN fluctuators. Thus we have achieved an exponential 
improvement in the efficiency of the noise approximation. 
Alternatively, we can choose the eigenvalues of T freely 
and obtain a valid matrix V with numerical optimization, 
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FIG. 1: Logarithm of the power spectral density for five 
independent RTN fluctuators (dash-dotted line), a multi-state 
Markovian source corresponding to 31 RTN fluctuators (solid 
line), and an ideal 1// noise (dotted line). The transition 
rates of the RTN fluctuators are in both cases distributed 
uniformly on the interval [70,3070]. 



which may result in even more faithful approximation of 
1// Q noise. 

Figure Q] compares the approximation of the spec- 
tral density of 1// noise generated by independent RTN 
sources and by a multi-state Markovian source. For the 
RTN approach, we choose 5 independent noise sources, 
for which the transition rates 7^ are uniformly distributed 
in the range [7min,7max] = [70,3070], and the strengths 
are given by = 1/^%. This yields a fiuctuator 
with 32 distinct noise states. For the multi-state fiuc- 
tuator, we choose a 32-state noise source, for which 
the nonzero eigenvalues Afe of its transition rate ma- 
trix r are distributed uniformly on [—6070,-270], and 
Xk = Afe/2. Thus the condition in Eq. ([9]) is satis- 

fied for both of the approaches and the multi-state noise 
source has an autocorrelation function and power spec- 
tral density which are equal to those for a certain ensem- 
ble of 31 RTN fluctuators. We employ representations of 
similar computational complexity here in order to be able 
to assess the relative accuracy for a given computational 
effort. 

Figure Q] shows that an ensemble of five RTN processes 
is not an accurate model for 1// noise, whereas a single 
32-state Markovian noise source is quite accurate, espe- 
cially in the range 370 < to < I670. The poor quality 
of the approximation with five RTN fluctuators is due to 
the small number of independent noise sources employed 
here, whereas the 32-state Markovian fiuctuator contains 
more parameters and thereby introduces more flexibility 
in the noise approximation. The frequency range over 
which the approximation is accurate is relatively short 
if one considers that the 1// noise detected in experi- 
mental applications often extends over several frequency 
decades. The width of this frequency range can of course 
be increased by increasing the width of the region from 
which the eigenvalues of the matrix T are chosen. In 
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this case, however, the number of discrete levels in the 
Markovian source must also be increased to preserve the 
desired accuracy. For the main purpose of demonstrating 
the feasibility of the numerical optimization algorithm, 
in the rest of this paper we will continue to approximate 
1// Q noise by a single Markovian noise source with 32 
levels. 



III. QUBIT DYNAMICS AND CONTROL 

In Ref . [24| , the temporal evolution of the qubit density 
matrix was calculated by averaging over 10 4 -10 5 unitary 
quantum trajectories, each corresponding to a sample 
noise path. To ensure accuracy, a large number of uni- 
tary trajectories are required, which results in extensive 
computational effort. In Ref. [29j, exact deterministic 
master equations describing the average temporal evo- 
lution of quantum systems under Markovian noise were 
derived. 

Following Ref. [29] , we introduce a conditional density 
operator pk(t) which corresponds to the density operator 
of the system averaged over all the noise sample paths 
occupying the fcth state at the time instant t. The con- 
ditional density operators are normalized such that the 
trace of the operator pk(t) yields the probability of the 
fcth noise state as -Pfe(t) = Tr [p k (t)}. The total average 
density operator can be expressed as 



(16) 



The dynamics of pk is obtained from the coupled master 
equations [29] 

d tPk (t) = 1 [H k (t),p k (t)]+J2^k j p j (t), (17) 



where H k (t) is the Hamiltonian of the system correspond- 
ing to the fcth noise state, and T k j the transition rate 
from the jth state to the fcth state, as defined in Sec. [ill 
Specifically, in our one-qubit case, 



H k {t) = \a{t)(j x + \b k a 



(18) 



where b k is the noise amplitude of the state fc. We shall 
use £ a {p} to denote the state p evolved under the influ- 
ence of noise and the control sequence a. 

The fidelity function quantifying the overlap between 
the desired state pi and the actual achieved final state is 
defined as 



<!>{pi,£a {po}) = Tr p\E a {p } 



(19) 



where po is the initial state of the system. To measure 
how close the evolution £ a is to the intended quantum 
gate operation U, we calculate the average of the fidelity 



4>(U poW ,£ a {po}) over all pure initial states po, and ob- 
tain the gate fidelity function 



(20) 



k—x,y,z 



We aim to find the optimal control pulses which max- 
imize the fidelity of the achieved quantum operation, 
and hence apply a typical gradient based optimization 
algorithm such as the gradient ascent pulse engineering 
(GRAPE) method developed in Ref. |311 |. If the contin- 
uous pulse profiles are approximated by piecewise con- 
stant functions, the gradient of the fidelity function with 
respect to these constant pulse values and durations can 
be calculated by the chain rule. This gradient is further 
used as a proportional adjustment to update the control 
pulse profile. The optimization procedure is terminated 
when certain desired accuracy is achieved. Note that 
due to the non-convex nature of the problem, the gra- 
dient based algorithm will only yield a locally optimal 
solution. We further employ a multitude of initial con- 
ditions to find a control pulse which achieves the highest 
fidelity. 



IV. QUANTUM MEMORY 

In this section, we focus on the implementation of 
quantum memory, i.e., the identity operator. For the pur- 
pose of comparison with the optimized pulse sequences, 
we introduce four other kinds of control schemes which 
generate the identity operator. 

The first reference sequence is simply not to apply any 
external control pulse, i.e., a(t) — 0. This pulse has 
no compensation for decoherence or error. The second 
reference sequence is a constant 2ir pulse given by 



for * e [0, 2irh/a n 



(21) 



The third reference sequence is the composite pulse se- 
quence known as compensation for off-resonance with a 
pulse sequence (CORPSE), which was originally designed 
to correct systematic errors in the implementation of one- 
qubit quantum operations and to provide high order con- 
trol protocols for systematic qubit bias, i.e., for the noise 
correlation time r c — » oo jH,!!^. For the identity oper- 
ation, the CORPSE pulse sequence can be obtained as 



Omax, for < t' < 7T 

asc27r(t) = { -aw, for ir < t' < 3tt 

fWx, for 37T < t' < 47T, 



(22) 



where the dimensionless time t' is defined as t' — a maK t/H. 

In the absence of noise, the CORPSE sequence gen- 
erates the identity operator exactly although it requires 
twice as long operation time as a 2ir pulse, the second 
reference pulse above. In the presence of small system- 
atic errors, the CORPSE sequence is much more accurate 
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than the 2-k pulse. For example, consider a state trans- 
formation from the north pole back to itself on the Bloch 
sphere. For rj(t) = A in Eq. (U), the fidelities defined in 
Eq. lfl9|) can be derived to be 

4 \ Umax / \ ^max / 

and 

0SC2. = 1 - 4tt 2 ( — ) +o(—) . (24) 

We observe that the error in the fidelity of the 2ir pulse 
is fourth order in the relative noise strength A/a max , 
whereas for the CORPSE pulse sequence it is eighth or- 
der. Thus the CORPSE sequence is much more accurate 
than a 2ir pulse in correcting the effects of systematic 
errors on quantum memory. 

The fourth standard pulse sequence which we take as a 
reference is the Carr-Purcell-Meiboom-Gill (CPMG) 0] 
sequence which is designed to preserve qubit coherence. 
In our context, this sequence consists of a ir/2 pulse fol- 
lowed by multiple ir pulses at intervals t p , followed by a 
final 7r/2 pulse to bring the system back to the original 
state. This pulse sequence is designed for T 2 measure- 
ments on spins, starting from the |0) state. Thus one 
does not expect a CPMG pulse sequence to perform as 
well if the initial state is averaged over the Bloch sphere 
as is done to compute a gate fidelity. 

We first present the fidelities obtained for the iden- 
tity operator using the various control pulse options in 
the presence of 1// noise. The noise is generated here 
by the single Markovian noise source discussed in Sec.ITu 
with transition rates distributed uniformly over the inter- 
val [l/r c , 30/r c ] . In Fig. O the fidelities obtained from 
optimized control pulses, 2-7T pulse, CORPSE, CPMG, 
and zero pulse sequences are plotted as functions of the 
characteristic correlation time r c of the approximate 1// 
noise. Here, CPMG1 and CPMG2 refer to two CPMG 
types of pulses with the intervals between ir pulses be- 
ing 7r and 2n, respectively. The total duration for these 
pulses are all 12irh/a max . The optimal control pulse is de- 
signed for 67r, and therefore we repeat it twice. Similarly, 
we repeat the 27r pulse 6 times, the CORPSE sequence 
3 times, the CPMG1 sequence 3 times, and the CPMG2 
sequence twice. The optimal control pulse yields clearly 
the highest fidelity among all these pulses, whereas the 
zero pulse sequence has the worst performance as there 
are no correction mechanisms. Note that due to motional 
narrowing, all curves approach unit fidelity in the limit 
t c ^0. 

The memory access rate is an important specification 
in modern computer technology [35j. In our context, it 
corresponds to the total duration of the control pulses. 
Figure [3] shows the fidelity as a function of the dura- 
tion for the numerically optimized control pulses. Equa- 
tion U) implies that in the absence of noise, the quantum 
system will generate an identity operator for a = a max 
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FIG. 2: Fidelity of the quantum memory as a function of 
the characteristic correlation time r c for optimized control 
pulses (black solid), a 2ir pulse (black dash-dotted), CORPSE 
pulse sequence (black dotted), CPMG1 pulse sequence (black 
dashed), CPMG2 pulse sequence (gray solid), zero pulse se- 
quence (gray dash-dotted). The operation time is chosen to be 
127rft/a max . The noise is produced by a single 32-state Marko- 
vian source with the average strength (\q\) = 0.125 x a max 
corresponding to 31 RTN fluctuators with the transition 
rates uniformly distributed over the region [l/r c ,30/r c ] and 
strengths chosen as described in Sec. HH 




T/(7r/i/a max ) 

FIG. 3: Fidelity of the quantum memory as a function of the 
operation time for control pulses optimized at each point. The 
noise is produced by a similar multi-state Markovian source 
as in Fig. [2] with r c = 3/i/amax. 



and the duration T = 2mr/a max . In Fig. 03 we observe 
that, despite an overall decrease, there are peaks in the 
fidelity near these operation times. Thus we can regard 
2mr/a max as the natural periods for quantum memory, 
and we always choose the total duration of control pulses 
correspondingly. 

Here, we study the relation between the optimized fi- 
delities achieved above and the average noise strength 
(\r)\), for a fixed value of the characteristic correlation 
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FIG. 4: Fidelity of the quantum memory as a function of the 
average absolute noise strength for optimized control pulses 
(black solid), 2tt pulse (black dash-dotted), CORPSE (black 
dotted), CPMG1 (black dashed), CPMG2 (gray solid), and 
zero (gray dash-dotted). The operation time is chosen to be 
127r/i/a max . Except for its strength, the noise is produced by 
a similar multi-state Markovian source as in Fig. [2] with the 
correlation time r c = 30/i/a max . 



time r c = 30ft/a max . Figure [4] shows the fidelity as a 
function of the noise strength for the optimized control 
pulses, 2tt pulse, the CORPSE, CPMG1, CPMG2, and 
zero pulse sequences. At small values of (\rj\) again, the 
optimized control pulses consistently achieve higher fi- 
delities than all reference pulses. However, we note that 
if the noise strength exceeds ^0.4, the optimized pulse 
sequence reduces to the zero pulse sequence, i.e., any 
nonzero pulse sequence will actually deteriorate the fi- 
delity performance. 

The discussion above is based on the specific noise den- 
sity spectrum l/f a with a = 1. Figure [5] shows the 
fidelities of quantum memory for four optimized control 
pulses, each of which is obtained for a different value of a. 
The noise is produced here by a single multi-state Marko- 
vian source with average strength (\rj\) = 0.125 x a max , 
and the total duration for all control pulses are fixed 
to 6ir. A systematic scaling of the correlation time axis 
with respect to a is clearly visible in Fig. [H This phe- 
nomenon is explained by the fact that the concentra- 
tion of the power spectrum of 1// Q to high frequencies, 
i.e., long correlation times, increases with a. Hence, the 
curves scale down in t c with increasing a. 



V. NOT GATE 

In this section, we focus on the generation of high- 
fidelity NOT gates, i.e., the u x operator, under 1// noise. 
As in the case of quantum memory, we compare the nu- 
merically optimized results with reference pulses. In this 
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FIG. 5: Fidelity of the quantum memory achieved with op- 
timized control pulses as a function of the characteristic cor- 
relation time t c for l/f a noise with a — 1 (solid ), a = 1.25 
(dotted), a = 1.5 (dash-dotted), and a — 1.75 (dashed). The 
operation time is chosen to be 67rft/a max . The noise is pro- 
duced by a similar multi-state Markovian source as in Fig. [2] 
with variable values of the power a. 

case, our first reference pulse is the it pulse given by 

for t € [0,irh/a max ], (25) 

which in the absence of noise is the most efficient way of 
achieving a NOT gate. In addition, we will use the two 
composite pulse sequences CORPSE and short CORPSE 
[321 . l33l | which assume here the form 

fflmax, for < t' < ir/3 
-amax, for 7t/3 < t' < 2tt (26) 
Omax, for 27T < t' < 137r/3, 

and 

{-a ma x, for < t' < tt/3 
amax, for tt/3 < t' < 2tt (27) 
-a max , for 2tt < t' < 7tt/3, 

respectively. Both of these pulse sequences correct for 
systematic error, CORPSE being more efficient. How- 
ever, the operation time of short CORPSE is much 
shorter than that of CORPSE, and hence it can yield 
higher fidelities in the presence of noise. 

Figure [6] shows the NOT gate fidelities obtained by 
the reference and optimized pulses in the presence of the 
same 1// noise as employed in the analysis of quantum 
memory in Sec. IIVI We observe that for long enough 
correlation times, the composite pulse sequences provide 
good error correction. Furthermore, as observed earlier 
for RTN [lij], for intermediate correlation times, short 
CORPSE achieves the highest fidelity among the refer- 
ence pulses. Figure [7] presents the pulse sequences ob- 
tained from the numerical optimizations for three dif- 
ferent values of the noise correlation time t c . For the 
optimized pulse sequence, we find a transition from an 
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approximately constant pulse to a short CORPSE -like 
pulse sequence at characteristic correlation time t c r* 
50fi/a max . This change in optimal pulse sequence is 
responsible for the apparent discontinuity in the first 
derivative of the fidelity curve in Fig. [6l 




200 



FIG. 6: NOT gate fidelities as functions of the characteris- 
tic noise correlation time r c for a it pulse (dotted), CORPSE 
(dash-dotted), short CORPSE (dashed), and gradient opti- 
mized pulse sequence (solid). The 1// noise is generated as 
in Fig. El 

These results for the generation of NOT gates under 
1// noise are qualitatively quite similar to the previous 
results presented in Refs. [24 l29| for a single RTN. This 
similarity is due to the fact that 1// noise can be re- 
garded as arising from a sum of independent RTN fluc- 
tuators, each of which having a similar fidelity depen- 
dence on their correlation times. Note that the scale for 
the reference correlation time t c of the fidelity obtained 
in presence of 1// noise in Fig. [6] is somewhat different 
from the corresponding scale for the correlation time of 
a single RTN source, since the 1// noise involves an en- 
semble of RTN fluctuators with a range of correlation 
times. 



VI. CONCLUSIONS 

We have studied a single qubit under the influence of 
l/f a noise for 2 > a > and investigated how deco- 
herence due to this noise can be suppressed in the im- 
plementation of single qubit operations. We presented 
an efficient way to approximate the noise with a discrete 
multi-state Markovian fluctuator. Due to this finding, 
the average temporal evolution of the qubit density ma- 
trix under 1/ f a noise can be efficiently determined from 
a deterministic master equation. 

Employing these exact deterministic master equations 
describing the temporal evolution of the qubit density 
operator under Markovian noise, we applied a gradient 
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FIG. 7: Optimized pulse sequences yielding the highest gate 
fidelities for correlation times (a) 45fi/a max , (b) 100?i/a max , 
and (c) 150fi/a max corresponding to Fig. [6] 



based optimization procedure to search for optimal con- 
trol pulses implementing quantum operations. In par- 
ticular, we studied the physical application of quantum 
memory, i.e., the identity operator, which is a funda- 
mental concept in the realization of a quantum com- 
puter. The optimized control pulses significantly im- 
proved the fidelity over several reference sequences such 
as 2tt, CORPSE, CPMG, and zero pulses. We observe 
peaks on fidelity curves corresponding to integer multi- 
ples of 2nh/a lrmyl in the total durations of control pulses, 
where a max is the maximum magnitude of the external 
control field. We also studied the performance of opti- 
mal control pulses under 1// Q noise for several different 
values of 2 > ev > 1 , and found a monotonic behavior in 
the noise frequency as a function of a, i.e., the fidelity 
curves are scaled down in the correlation time for increas- 
ing a. We also investigated how the fidelities degraded as 
the noise strength increases. For the generation of high- 
fidelity NOT gates, we obtained results showing qualita- 
tively similar behavior to the previous results presented 
in Refs. [H, H§] for a single RTN source. In particu- 
lar, just as for a single noise source, in the presence of 
1// Q noise we observed a transition in the optimal control 
pulse sequence from a constant pulse to a CORPSE-like 
sequence as the noise characteristic correlation time r c is 
increased. 

This approach of coupled master equations indexed by 
noise states of the environment, together with an opti- 
mization technique for pulse design can be readily gen- 
eralized to multiple qubits evolving in the presence of 
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1// Q noise and other Markovian noise sources. Further- 
more, it can be used to develop realistic pulse sequences 
for mitigation of nuclear spin and surface magnetic noise 
acting on donor spins implanted in silicon [111], as well 
as for suppression of background charge noise acting on 
superconducting qubits Q. In future, we will study the 
implementation of multi-qubit gates, e.g., the controlled 
NOT gate, in noisy systems and the swapping of quan- 
tum information from a noisy qubit to long term quan- 
tum memory. We will also consider more realistic noise 
with l/f a spectrum over many frequency decades. 
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